% 得到最大的非0基频
function baseFreq = getMaxBaseFreq(x, fs)
    n = length(x);
    k = 0:floor(n/2);
    y = fft(x);
    f = k * fs / n;
    p1 = abs(y(k+1)) * 2 / n;
    f_nonzero = f(2:end);    % 非0频率轴
    p_nonzero = p1(2:end);   % 非0频率对应的幅度
    [~, max_idx] = max(p_nonzero);  % 最大幅度的索引
    baseFreq = f_nonzero(max_idx);  % 对应的非0频率
end